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Nuclear effects in CC interactions with carbon nuclei have been studied by using backward going protons 
and n~ . Detailed analyses, of the momentum distributions and of the production rates, have been carried out in 
order to understand the mechanism producing these particles. The backward proton data have been compared 
with the predictions of the reinteraction and the short range correlation models. 



1. Introduction 

Since a long time p| it has been observed that 
in the high energy interactions off nuclei there 
are particles emitted backwards, with respect 
to the beam direction, which have energies not 
allowed by the kinematics of collisions on a free 
and stationary nucleoli. Backward going pro- 
tons are commonly observed while, in absence 
of nuclear effects, their production is forbidden. 
High energy mesons, whose production in the 
backward direction is only allowed up to a given 
momentum, are detected also, at momenta above 
such a limit. Essentially two types of models 
have been proposed to explain the origin of these 
particles. The first is based on the intranuclear 
cascade (INC) mechanism. The production of 
particles in the kinematically forbidden region is 
seen as the result of multiple scattering and of 
interactions of secondary hadrons, produced in 
the primary z^-nucleon collision, with the other 
nuclcons while they propagate through the nu- 
cleus |U^]. The second approach explains the 
production of backward particles as the result 
of collisions off structures with mass larger than 
the nucleon one. These structures are formed, at 
small interparticle distance, under the action of 
the short range part of the nuclear force. They 
may be described as clusters of a few correlated 
nucleons pj] or multiquarks clusters ||[3j. The 
nucleons in these structures can acquire high 
momenta and the fast backward going particles 
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can be seen as a direct manifestation of the high 
momentum tail of the Fermi distribution. 
These two classes of models are not mutually 
exclusive and the production in the backward 
hemisphere can have contributions from both 
mechanisms. 



2. The NOMAD detector and event selec- 
tion 

The NOMAD detector § consisted of an ac- 
tive target of 44 drift chambers, 3 x 3m 2 each, 
of low average density (0.1 g/cw?) and overall 
mass of 2.7 tons. The drift chamber composi- 
tion was made by over 90% of the total weight by 
carbon and elements with nearby atomic num- 
ber. Downstream the beam the target was fol- 
lowed by a transition radiation detector which 
provided electron identification, a preshower and 
an electromagnetic calorimeter. All these systems 
were located inside a dipole magnet which created 
a magnetic field of 0.4 T. A hadron calorime- 
ter and two muon stations were located just be- 
hind the magnet coil. NOMAD collected data 
from 1995 to 1998. During the run a total of 
5.1 x 10 19 protons on target, corresponding to 
about 1.3 x 10 6 ^, CC interactions, were recored 
on tape. The events are selected by requiring one 
primary vertex with at least two tracks, one of 
which has to be a negative muon of momentum 
greater than 3 GeV/c. The sample selected for 
this study amounts to 944019 events. 
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Figure 1. Distributions of length vs. momentum for positive (left) and negative (right) backward going 
tracks selected as described in the text. The line indicates the position of the cut. 



3. Backward particle identification 

3.1. Backward proton identification 

The separation of the backward going pro- 
tons (Bp) from the other backward going posi- 
tive tracks (essentially n + ) is done by exploiting 
the range-momentum relation. For this reason 
the candidate tracks are required to range out 
in the drift chamber volume. Fig. [j] shows the 
experimental distributions of length vs. momen- 
tum for positive (left) and negative (right) tracks, 
going backward with respect to the beam direc- 
tion, and satisfying our selection criteria. Two 
distinct populations are clearly visible on the pos- 
itive sample plot. Protons, having a shorter range 
than 7r + , tend to accumulate in the lower right 
part of the plot while the 7r + 's tend to popu- 
late the left-hand side. Comparing the two plots 
we see that the lower right part of the negative 
sample is much less populated than the corre- 
sponding region of the positive one, these tracks 
being mainly n~ . Any positive backward going 
track with length smaller than the cut (the line 
on Fig. is identified as a Bp. The identified 
sample is then corrected for reconstruction ef- 
ficiencies and for tt + contamination (s=s 8% for 
Prec > 250MeV/c). For further details on the 
identification and correction procedure refer to 



3.2. Backward 7r identification 

In the case of backward going tt~ {Bit~) it is 
not necessary to look for tracks ranging out in 
the target. The choice of the negative charge 
allows for a sample of high purity; it is shown 
by MC studies that the e~ contamination is at 
the level of few percent above a momentum of 
« 200 MeV/c. We therefore assume that any 
negatively, backward going, charged track with 
Prec > 0.2 GeV/c is a n~ . The identified sam- 
ple is then corrected for reconstruction efficiency. 

4. Backward p and 7r invariant momen- 
tum distributions 

The inclusive spectrum of backward particles is 
typically represented using the normalized invari- 
ant cross section (1 /error) E d 3 a/dP 3 , where E 
is the energy of the backward going particle. The 
invariant cross section is usually parametrized by 
an exponential form as: 

1 E dN 
P dp* 



C e 



-BP* 



(1) 



where N ev is the total number of events, C is a 
constant and B is the slope parameter. Early 
hadronic experiments had found that B is al- 
most independent of the projectile type and mo- 
mentum, and of the atomic number of the tar- 
get. This behaviour has been termed "nuclear 
scaling". Neutrino experiments p|-p3[ have con- 
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Figure 2. Invariant momentum distributions for backward going protons (a), ir~ (b) in v^CC events and 
protons in v^CC events (c) 



firmed these properties. 

The inclusive spectrum of Bp and Btt~ is shown 
in Fig. and |^b respectively, together with the 
exponential fit. In the Btt~ case we have not in- 
cluded the first two points in the fit. These points 
are in a momentum region where, backward pro- 
duction being kinematically allowed, there are ad- 
ditional contributions not coming from nuclear ef- 
fects. The values of the measured slope param- 
eter B are reported in Table The first and 
the second errors are statistical and systematical, 
respectively. The systematic uncertainty was es- 
timated by slightly changing the values of all the 
cuts used, by varying by a small amount the cor- 
rection functions and also by changing the size of 
the fiducial volume. The values of the slope pa- 
rameters measured in this experiment are found 
to be compatible with the results obtained in 
other neutrino and hadronic experiments. The 
invariant cross section for Bp is larger than the 
one for Btt~ by about one order of magnitude but 
the values of the slopes are similar. The kinematic 



ranges of the two invariant distributions are also 
different. To be identified Bp have to stop in- 
side the target volume; the rather low density of 
the NOMAD target restricts to w 0.7 GeV/c the 
maximum useful momentum value. 
The invariant cross section and slope parameter 
for Bp in CC events are given in Fig. ||c and 
Table |l| respectively. During normal operations 
CC events were also collected due to the small 
Pfi component of the dominantly beam. A 
dedicated v run yielded an additional sample of 
Dft CC events included for this analysis. 
Antineutrino events are selected under the as- 
sumptions that the efficiencies and pion contami- 
nation are the same as those used for the neutrino 
events, but requiring a positively charged muon 
instead of a negative one. The final CC sample 
consists of 61 1 34 events. 
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Table 1 

Fit ranges and values of the slope parameter B, for backward protons and tt~ , as obtained from the expo- 
nential fit to the invariant cross section. The first and the second errors are statistical and systematical, 
respectively. 







AP (GeV/c) 


B(c 2 /GeV 2 ) 


C(c 3 /GeV 2 ) 


Backward p 


v M CC 


0.37 -=-0.70 


10.54 ±0.20 ±0.5 


4.08 ±0.19 ±0.5 




^ cc 


0.37 -=-0.70 


10.79 ±0.78 


2.71 ±0.54 


Backward 7r~ 


v„ CC 


0.32 -=-0.85 


10.03 ±0.28 ±0.3 


0.17 ±0.01 ±0.02 



5. Energy dependence of the slope param- 
eter 

In Fig. |3| we show the slope parameter B for 
protons as a function of the visible hadronic en- 
ergy Eh ad and the Q 2 of the event. The visible 
hadronic energy is defined as: 

Eh ad = ^ E c + 2^ En 

where Yl E c is the sum of the energies of re- 
constructed charged tracks (assuming the mass 
of the pion if the particle is not otherwise iden- 
tified) and E n includes tracks associated to 
secondary vertices corresponding to interactions 
of neutral particles and the energy of neutral 
particles reconstructed in the ECAL. Q 2 is the 
square of the four-momentum transfer, Q 2 = 
4 E vis E ll sin 2 8/2, where 9 is the muon angle with 
respect to the neutrino direction. No significant 
dependence of the slope B on either quantity is 
observed, in agreement with the expectations of 
"nuclear scaling" as observed in hadronic beam 
experiments. The range covered by NOMAD is 
similar to the one covered by different experi- 
ments with hadronic beams. 



6. Backward particle rates 

In order to study a possible atomic number de- 
pendence the Bp rate has been compared with 
the results of the other i/-nucleus experiments. 
The measured yields have all been normalized to 
the same momentum interval (350 to 800 MeV/c) 
assuming the dependence given in Eq. [j] with the 
measured slopes. 
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Figure 3. The value of the slope parameter B 
shown as a function of the hadronic energy Eh ad 
(top) and of Q 2 (bottom). 



The extrapolation for NOMAD gives: 

< Nb p >350-H800AfeV/c = 

[ 52.8 ± 0.6{stat.) ± 7(syst.) ] x 10~ 3 

The A dependence for experiments most directly 
comparable to NOMAD is shown in Fig. [|a. In 
the range A = 20 -=- 80 it has been parametrized 
as < N B p >oc A a , where a = 0.68 ± 0.12 JT|. 
The same parametrization with a similar power 
law was found to describe Bp production in tt + 
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7. Comparison of data with models 
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Figure 4. The average number of Bp per event (a) 
and of Bir~ per event (b) in the momentum range 
from 350 -j- 800MeV/c in neutrino experiments 
as a function of the atomic number A. The line 
shown in (b) is the result of the fit described in 
the text. 



7.1. Eh ad and Q 2 dependence of Bp and 
Bir~ rates 

The Bp and the Bir~ rates have been stud- 
ied as a function of the hadronic energy Eh ad 
and of Q 2 . In both cases, shown in Fig. ||, a de- 
crease of the yield with increasing Eh ad and Q 2 
is observed. This can be interpreted in terms of 
the "formation zone" concept. This is the dis- 
tance (or the time) from the production point 
which is required for the secondary hadrons to 
be "formed" , i.e. to be able to interact as phys- 
ical hadronic states. Since the distance/time re- 
quired, due to the Lorentz time-dilation factor, is 
proportional to the energy of the secondary, the 
INC process is restricted to slow hadrons which 
have formation lengths smaller than the nuclear 
radius. The larger Eh ad and Q 2 , the larger is 
the average energy of the outgoing partons there- 
fore resulting in hadrons which have higher prob- 
ability to be formed outside the nucleus. As a 
consequence, reinteractions will decrease and so 
will the slow proton rates. In the Btt~ case the 
dependence of the yield on Q 2 and on Eh ad is 
less pronounced. This can be understood since 
the Btt~ rates are a less sensitive probe of nu- 
clear effects because Bir~ production on a sta- 
tionary nucleon is kinematically allowed for mo- 
menta up to about half the nucleon mass. Fur- 
thermore Btt~ can be produced in the decay of 
forward going resonances or unstable particles. 



and K + collisions with Al and Au nuclei at 250 
GeV/c [l4 |. It is evident from Fig. ||a that this 
simple power law does not describe the NOMAD 
data taken at a lower < A >. The Btt~ rate was 
directly measured in the same momentum range. 
Its average value is found to be: 



Btt- >350^800AfeV/c 



-3 



[ 4.8 ± 0.l(stat.) ± 0.3(syst.)] x 10 



In this case a good fit to the two BEBC points 
and to the NOMAD value can be obtained using 
the form A a giving a — 0.83 ± 0.25, as shown in 
Fig.@b. 



7.2. Effects of short range correlations in 
Bp production 

According to the picture proposed by the cor- 
related nucleon/quark models the Bp production 
is explained as the result of a neutrino interaction 
within a correlated cluster of nuclcons/ quarks. 
This cluster is formed, for a short time, when two 
nucleons in their motion inside the nucleus ap- 
proach each other so closely as to come under the 
effect of the short range component of the nuclear 
force (r c = 0.5 -r- 0.7/m). The repulsive charac- 
ter of this interaction is responsible of the high 
relative momenta of the cluster members. The 
simplest, and the most probable case when con- 
sidering the overlapping probabilities of nucleons 
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Figure 5. Top picture: the average number of Bp (370 < P < 700MeV/c) per event as a function of the 
hadronic energy Eh ad fleft) and of Q 2 (right). Bottom picture: the average number of Bir~ (350 < P < 
800MeV/c) per event as a function of the hadronic energy Eh ad (left) and of Q 2 (right). 



in the nucleus, is the two-nucleon case. Accord- 
ing to this model the Bp is released when the 
correlated pair is broken in the ^-cluster interac- 
tion. If the effects of reinteractions are neglected, 
the released Bp can leave the nucleus keeping its 
original momentum. To study these correlations 
it is customary to use the variable a defined as: 



[E -Pi)/M 



(2) 



where E and Pi are respectively the energy and 
longitudinal momentum of the Bp and M is the 
nucleon mass. For Bp, a > 1 since Pi is negative. 
In this model, due to the target motion a correla- 
tion between a and the variable v = xy, where x 
is the Bjorken scaling variable x = Q 2 /2MEhad 
and y = Ehad/(Ehad + E^) is expected. The 



variable v is related to the muon kinematics and 
can be written as: 



(E„ - Pi)/M 



(3) 



where E^ and P 1 ^ are the muon energy and longi- 
tudinal momentum. The (v, a) correlations were 
searched for in the data by calculating for each 
event a and v as defined in Eq. and |^. For each 
a bin we plot the variable < vn > defined as: 



<v N > = 



< V > B p 



< V > no Bp 



(4) 



where < v > no Bp is obtained from the full sam- 
ple of events without a Bp. According to Ref. Q 
the average values of v in events with a Bp, 
< v >Bp, is related to the average value of v in 
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events where no Bp is present, < v > no Bp-, by: 

< V > B p = <V > no Bp (2 - a) (5) 

More generally for a cluster composed of £ nucle- 
ons the relation is ||: 

f a\ £ 

< V > B p = <V > no Bp I 1 - 7 I |— ■ - (6) 

Fig. ^a shows < % > as a function of a. The 
data indicate a slope of —0.22 ± 0.06 with a 
X 2 /ndf = 5.3/6 (the x 2 /ndf in the hypothesis of 
no dependence on a is 19.1/7). The two-nucleon 
correlation mechanism fails to describe our data. 
Either higher order structures are playing a lead- 
ing role H or the observed low level of correlation 
is due to the presence of reinteraction processes. 
In this case part of the Bp are emitted as a result 




Figure 6. The variable < Vjsr > plotted as a func- 
tion of a. The lines represent the predicted cor- 
relation (Eq. ^) for a number £ of nucleons in 
the cluster equal to 2, 3 and 4- In (a) all the 
Bp events were used; in (b) only events having 
a Bp with cosOj < 0. 6j being the angle of the 
Bp with respect to the hadronic jet direction. 

of reinteractions in the nucleus and are not related 
to the target nucleon. The presence of intranu- 
clear cascade processes could therefore dilute the 
existing correlation to the observed level. To test 
this hypothesis we tried to reduce the component 
due to rescattering in the selected Bp sample. 
Having observed the correlation existing between 
the multiplicities of slow tracks and rescattered 



protons we applied increasingly tighter cuts on 
the number of slow tracks (P < 700 MeV/c). As 
a consequence of these cuts the degree of correla- 
tion between a and v increases. The fit values of 
the slopes are reported in Table || together with 
the definitions of the cuts applied. We have also 
observed a strong correlation between the pres- 
ence of protons travelling backward in the lab but 
forward with respect to the hadronic jet direction, 
and the concentration of events at small Q 2 val- 
ues and large angles with respect to the beam. 
Since also a small Q 2 indicates the presence of 
rescattering, the exclusion of these events should 
highlight the expected correlation. The result- 
ing slope is -0.39 ± 0.07 with a x 2 jndf = 4 -9/6 
(see Fig. ^|b). The observed behaviour is consis- 
tent with the hypothesis of the correlations effects 
being to some degree hidden by the presence of 
rescattering. Reducing the rescattering compo- 
nent these correlations seem to become stronger. 

Table 2 

The fitted value of the (a,v) slope, and the cor- 
responding x 2 l n df , for Bp selected from events 
with various numbers of positive (n + ) and nega- 
tive (n~ ) low momentum (P < 700MeV/c) par- 
ticles. 



n + n 


slope 


X z /ndf 


< 3 < 2 


-0.23 ±0.06 


5.7/6 


< 2 < 1 


-0.25 ±0.07 


3.3/6 


< 2 


-0.30 ±0.07 


3.7/6 


1 


-0.37 ±0.10 


3.9/6 



8. Comparison with the NOMAD Monte 
Carlo 

The NOMAD event generator NEG-N is based 
on a modified version of LEPTO 6.1 and 
JETSET 7.4 fl6| (with fragmentation parame- 
ters tuned on the NOMAD data) for the deep 
inelastic simulation and on dedicated generators 
for quasielastic events and events in the resonance 
region. Nuclear effects are taken into account 
by incorporating the formation zone intranuclear 
cascade (FZIC) code of DPMJET-II.4 @. For 
the hadrons produced inside the nucleus the code 
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Figure 7. Invariant momentum distributions of 
Bp simulated with NEG-N for different values of 
the formation time parameter. 

performs a complete simulation of tracking and 
of interactions in nuclear matter developing the 
cascade through several generations. The Fermi 
momentum is simulated according to p8|j . The 
event generator used for the oscillation search did 
not take into account the INC. In that case a dis- 
crepancy of 200MeV/ c on the average transverse 
missing momentum, p™ (the difference between 
the reconstructed hadron jet and muon transverse 
momenta), was found between the data and the 
MC prediction. In view of its importance in the r 
search, based on the exploitation of the different 
kinematics of the v r interactions with respect to 
the v e and ones in the transverse plane, the 
MC sample could not be directly used to calcu- 
late background and signal efficiencies. However, 
a method [jl9| was succesfully developed in or- 
der to correct the MC predictions using the data 
themselves. The results of this upgraded gener- 
ator show how the INC degradates the hadronic 
energy by producing n and slow p (emitted at 
large angles and difficult to detect). The simula- 
tion of nuclear effects reduces by a factor of two 
the p™ discrepancy between data and MC; the re- 
maining effect is probably due to a too optimistic 
detector simulation. In the code the hadron for- 
mation length is left as a free parameter to be 
supplied by the user. It has been tuned on the 



NOMAD data in order to achieve the best agree- 
ment between the simulated and the experimen- 
tal hadronic distributions of multiplicities, mo- 
mentum spectra and angular distributions. The 
best agreement is found for a formation time of 
2 fm/c. For this value the MC predicts a rate 
of < N Bp >350+800Mevy c = 53.0 x 10~ 3 with a 
slope of 10 in excellent agreement with the exper- 
imental results. It was found also (see Fig. |?]) that 
the rate is strongly dependent on the formation 
time while the slope is not. 
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